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Abstract 

Large long-distance standard model effects in flavor-changing neutral current (FCNC) 
semileptonic D decays can make observable these processes in future measurements. 
Eventual disagreements in this sector and/or the observation of lepton family violating 
(LFV) D decays would require an explanation beyond the standard model framework. 
In this paper we confront present experimental data on leptonic and semileptonic FCNC 
and LFV D meson decays with a version of the two-Higgs doublet model that allows 
these effects to occur at tree-level. The stringent bounds on the parameters of the 
model are obtained from — > l'^l'~ and D — > nl^l'^ decays. The consistency of the 
model requires that the branching fractions oi D ^ Vl^l'^ decays should be below the 
10-9 level. 



PACS: 12.15.Mm, 12.60.Fr, 13.20.Fc 



1. Introduction 

Flavor changing neutral currents (FCNC) in leptonic and semileptonic decays of charmed 
mesons are, higher order, very suppressed modes in the standard model (SM) of particle inter- 
actions mi-p. The short-distance contributions to these processes in the SM are expected to 
give branching fractions at the 10~^^ level for —>■ fj^^fi^ and 10"^ for D nl^l^ processes, 
while long-distance effects can enhance these predictions up to 10^^^ |^ and 10^'' ~ 10^^ 0, 
respectively. Present experimental upper limits for these decays are in the range 10^^ ~ 10"^ 
for l+l'- and lO^^ ~ lO^^ for D Xl+l'- 0-0 (X is a pseudoscalar or vector 

meson and /, /' = e, yu). On the other hand, lepton family violating (LEV) processes, i.e. 
/ 7^ are completely forbidden in the SM scenario with unmixed lepton generations. Thus, 
FCNC and/or LEV leptonic and semileptonic D decays can serve to test the mechanisms 
responsible of long-distance contributions or eventually would require an explanation beyond 
the SM framework. Yet another (unlikely) possibility is that nature places FCNC processes 
well below the SM expectations. This would force to revise the estimates of long-distance 
effects or, again, to invoke beyond the SM contributions to explain the eventual destructive 
interference with the SM amplitudes. 

Recently, the study of FCNC in charm quark decays has attracted a renewed interest 
m^llll- hand, it has been pointed out that these rare decays in models of new 



physics can be enhanced over the SM predictions by several orders of magnitude [Q. On the 
other hand, the existing bounds on FCNC and LFV D decays have been improved recently at 
FERMILAB E791, E771 and E687 experiments and by the CLEO Collaboration [ig . 

In addition, some projects have been proposed with the aim to reconstruct of the order of 
10^ charm decays during the Tevatron Run II |TT]| , which would increase the sensitivity 
to FCNC and LFV processes by almost three orders of magnitude with respect to present 
experiments. Therefore, it becomes timely to explore all possible scenarios of new physics 
that may give sizable contributions to these rare decays. 

In this paper we consider the constraints imposed by FCNC and LFV D meson decays on 
a general two-Higgs doublet model that allows these effects to contribute at tree level . 
The variant of the model considered here is built in such a way that tree-level FCNC inter- 
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actions of the neutral Higgses do not spoil the good agreement between the SM predictions 
and experiment for the down quark sector. The constraints on Yukawa interactions of the 
charged Higgses of this model have been studied in previous works [|T^. Here we consider 



the effects of Yukawa interactions of the neutral Higgses in FCNC and LFV decays of D 
mesons. To be more specific, we study the effects of neutral Higgses of this model in the 
D° l^l'^, D Pl^l'^ and D Vl^l'^ decays (-P(V^) stands for a pseudoscalar (vector) 
light meson and 1,1' = e or /i), which will provide a rather wide set of constraints on the 
effective Yukawa couplings of the model. 

2. The model. 

The variant of the two-Higgs doublet model needed in our work has been described 
elsewhere |]13[. The general form of the Yukawa interactions that allows tree- level FCNC 



processes is given by |12 



+'¥[{K^2 + ^K'^i)ll + h.c., (1) 

where F, F', G, G',K and K' are dimensionless 3x3 matrices, Ql = {Ul, Dl) with f/° 
{DD the triplet of left-handed up (down) quarks, and \E'° = {^l^ ^l) has a similar definition 
in terms of leptonic fields. C, parametrizes the small breaking of the discrete symmetry that 
forbids FCNC at tree-level. The superscript in fermion fields stands for weak eigenstates. 

Since we are interested in having FCNC contributions only in the up-quark sector, we 
shall drop the term proportional to G' in Eq. (1) [|13|. Notice that the Yukawa interactions 
for leptons are built to allow FCNC in the charged leptons and keep massless neutrinos. After 
spontaneous symmetry breaking, with = {0,vi/\/2) and ($2)"^ = (0, f2e~*°'/v^), the 

model contains five physical Higgses; the mass matrices for quarks and charged leptons 
become: 

Mu = ■^{Fv, + ^F'v2e-'^') , (2) 
Md = ■^Gv2, (3) 
Ml = l={Kv2 + ^K'vie-'"') . (4) 



For simplicity we choose to work in a basis where Mu and Mi are diagonal. Notice that, 
unlike the case where ^ = 0, F and F' (respectively, K and K') are not diagonal matrices 
and can allow for (unsuppressed by fermion masses) FCNC interactions in the up-quark 
sector. 

The Yukawa interactions between mass eigenstates of neutral scalar Higgses {Hq and ho), 
the pseudoscalar Higgs (Aq) and the fermions {U — {u, c, t) and I — {e,ii,T)) are given by 
(we do not write the interactions of down quarks because we are interested in FCNC in the 
up sector): 

Cn — ^t/{(FcosQ; + ^F'sinQ;)i7o + (— -P'sino; + ^F'cosQ;)/io 
v2 

+i{F sin /3 - ^F' cos /3) A075} U 

H — {{Ksma + ^K' cosa)HQ + {Kcosa — $,K' sma)ho 
v2 

+i{K cos (3 - iK' sin /3)^o75} I ■ (5) 

In these expressions, a is the angle that appears in the diagonalization of the neutral 
scalar Higgs bosons and tan/? = V2/vi . 

Due to the low energy scales involved in charm meson decays it becomes convenient to 
write out an effective four- fermion interaction Hamiltonian to describe the tree-level processes 
of our interest. The form of this Hamiltonian is: 

Heff = ^ {uAh.U ■ ILhJ + UAhoU ■ iLhJ + UAao^'U ■ lL^o75^} ■ (6) 

Using Eqs. (2) and (4), the effective couplings Aj and Li can be written as: 

^ ^^^F'(sinQ;-tan/?cosQ;e-^"') , (7) 
grriHo 

^ '^mw^^,.^ a + tan /? sin ae"'"') , (8) 
grriho 

^ _^!!!l^l^^'(cos/3 + tan/3sin/3e-'"') , (9) 
grriAo 

\/2Misma '^f^w .a ■ -ia'\ nn\ 

Lffo — n -\ 4 A (cos a — cot p sm ae ), (10) 

niHo sm/? gniHo 

a/2M/cosq; 2mH/^„,. . , ^ _,v\ /-.-.x 

Lho = ^ -^K'(sina + cot/3cosae , (11) 

uiho sm(3 grriho 

rriAo gruAo 
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As already anticipated, the leptonic couplings contain a (diagonal) piece proportional 
to fermion masses^ and another (non-diagonal) piece which is not a priori suppressed by 
fermion masses and will induce FCNC interactions. 

If we assume a specific ansatz for the Yukawa couplings F' and K' , we can use the exper- 
imental data on D decays to get bounds on the remaining parameters of the model. Instead, 
in the following we choose to use the available data to constrain the effective couplings given 
inEqs. (6)-(12). 



3. Constraints from leptonic and semileptonic D decays. 

The relevant hadronic matrix elements of the uc and M75C currents can be computed 
from the divergence of the c ^ d vector and axial vector charged currents and using isospin 
symmetry. Thus, we obtain: 

{Q\u^,c\D\p)) = ifj, , (13) 

(7r+(p')|wc|D+(p)) = V2{t:\p')\uc\D\p)) , (14) 

- f!^i^)Fr-(.^), (15) 

(y{p\e*)\uc\D{p)) = 0, (16) 

(p+(p',e*)|^75c|D+(p)) = V2{p\p\e*)\u^,c\D\p)) , (17) 

^'""^ -q.e*Af--'{q^) , (18) 



where q = p — p' is the momentum transfer to the lepton pair and e* is the polarization 
four-vector of the outgoing vector meson. In Eq. (16) \^ is a vector meson. Notice that the 
matrix elements for the D P and D ^ V transitions depend on only one form factor at 
the time. This happens because only the relative wave / = and / = 1 of the P-Higgs and 
\^-Higgs systems contribute to these transitions, respectively. 

For the D meson decay constant we take the value fo = 217 MeV which is obtained 
from the relation fn/ fos ~ 0.9 [|1^ and = 241 MeV from |T^. The g^-dependence of the 
scalar and pseudoscalar form factors appearing in Eqs. (14)-(18) are chosen to be monopolar 

Foiq') = , ^ Mq') = , (19) 

1 — g^/mp+ l — q-'/m^- 



^ Since we are interested in c — > m transitions we do not write a corresponding diagonal mass term in the 
quark couplings. 
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where mo+ and mg- are the masses of the scalar and pseudoscalar neutral D mesons, re- 
spectively. The normalizations of these form factors at = are taken from the relativistic 
quark model of Wirbel-Stech-Bauer ||16|| . 



The other hadronic matrix elements needed in our calculation are fixed either by identi- 
fying the uu content of final state isosinglet mesons, namely: 

{rj\uc\D^) = -^{7r°\uc\D°){cosep-V2smep), (20) 
v3 

{u;\uj5c\D'>) = (pV75c|/^°) , (21) 

or using SU(3) flavor symmetry: 

{K-^\uc\Dj) = {tt+\uc\D+) , (22) 
{K*+\u-f,c\Dj) = {p+\u-f5c\D+) . (23) 

Notice that we assume ideal uo — cf) mixing and we use 9p = —20° in Eqs. (20)-(21). 

The information on the experimental data about the FCNC and LFV D decays is taken 
from the 1997 update of Ref. Q, which already incorporates some recent results of Refs. 



-m 



In Table 1 we show the upper bounds for the products of couplings constants that can 
be constrained from the experimental data considered. We have introduced in Table 1 a 
short notation for coupling constants. First, we express the bounds from leptonic decays 
and D Vl^l'^ decays in terms of a''' = A^'j^^L^'^. Since both scalar neutral Higgses 
contribute to D ^ Pl^l'^ we have expressed the upper bounds in terms of the quantity 

Despite the fact that all the upper limits on branching ratios are at the 10~^ ~ 10~^ level, 
the different bounds on the effective couplings spread over two orders of magnitude. From 
Table 1 we conclude that the stronger bounds on the a''' couplings come from purely leptonic 

decays, while the same bounds from D Vl^l'~ decays are rather weak. Therefore, 
in the context of the present model, the leptonic decays imply that branching ratios of 
three-body decays of D's involving vector mesons should be below the 10~^ level. On the 
other hand, the best constraints on the a"' couplings are obtained from the D -kI^I'~ 
mainly because of the phase space suppression in the decays involving the rj meson. Finally, 
since the y-Higgs system in D ^ Vl^l'^ decays is in a / = 1 relative wave, this gives a 
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further phase space suppression and the absolute numerical bounds on the a'^ 's becomes 
weaker than the limits on the cr''''s (obtained from D P transitions). 

In order to draw any information on the Yukawa couplings of our interest let us make 
some considerations. To start, let us neglect the first term^in Eqs. (10)-(12) and set a' = 0. 
In this case we obtain the following expressions for a"' and a''': 

v^G^m^g sin /5 cos /5 
^ 1 i^FTi^KY i sm\a - f3) ^ cos\a - (3) \ ^ 

-J2Gf sin /5 cos /3 \ J 

or the relationship 

In the absence of information regarding the parameters of this model we will assume 
tan/3 ~ 1, rriho — ^30 GeV and rriHo — = 300 GeV. From Eq. (24) and the bounds on 
a'^' obtained from leptonic decays (see Table 1) we derive: 

(^^O'^l^^')"" < 2.9 X 10"^ (27) 
i^F')''%^KY^' < 1.7 X 10-3 (28) 
(^^')"^(^^')^' < 3.6 X 10^1 (29) 

Therefore, one may conclude that present experimental data on FCNC and LEV D decays 
only mildly constrain the strength of products of the relevant Yukawa couplings of this model. 
Since the (diagonal) terms proportional to fermion masses in Eqs. (10)-(12) are of 0(10—^) 
for the D Xfx'^fi' modes, the approximation done to derive Eqs. (24)-(25) is justified in 
view of the present experimental upper limits. 



Note that if a specific ansatz is assumed for these Yukawa couplings [|T7|, then Eq. (24) 
can furnish the allowed region for function of f3. Let us notice however, that Eq. 

(25) does not provide additional constraints on the Yukawa couplings unless, in addition, 
some information on the mixing angle a is introduced by hand. 

In summary, in this work we have studied the constraints imposed by FCNC and LEV 
leptonic and semileptonic D decays on a version of the two-Higgs doublet model that con- 
tains these effects at tree-level. The stringent bounds on the relevant Yukawa couplings are 
^Notice that this approximation is not necessary in the case of LFV decays. 
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obtained from two-body leptonic D° decays which are mediated by the pseudoscalar Higgs 
boson of the model. The best constraints on the Higgs scalar interactions are obtained from 

D 7rl~^l'~ decays. The three-body D decays involving vector mesons provide only very 
weak bounds and their measurements would have to be improved by five orders of magnitude 
in order to furnish similar constraints on the model as obtained from purely leptonic decays. 

Acknowledgements. The authors are grateful to Conacyt (G.L.C.) and Colciencias 
(R.M. and J.H.M.) for financial support. 



7 



References 

[1] K. S. Babu, X-G. He, X-Q. Li and S. Pakvasa, Phys. Lett. B205, 540 (1988) 

[2] A. J. Schwartz, Mod. Phys. Lett. A8, 967 (1993). 

[3] G. Eilam, J. L. Hewett and A. Soni, Phys. Rev. D44, 1473 (1991). 

[4] S. Pakvasa, e-print |hep-ph/9705397| and references cited therein. 



[5] P. Singer and D.-X. Zhang, Phys. Rev. D55, R1127 (1997). 

[6] Particle Data Group, R. M. Barnett et a/., Phys. Rev. D54 (Part I) (1996). See also the 
1997 update available at [http : //pdg . Ibl . gov A - 



[7] E791 Collaboration, E. M. Aitala et al, Phys. Rev. Lett. 76, 364 (1996). 

E771 Collaboration, T. Alexopoulos et al, Phys. Rev. Lett. 77, 2380 (1996). 
[9] E687 Collaboration, P. L. Fabretti et al, Phys. Lett. B398, 239 (1997). 



[10 

[11 
[12 

[13 

[14 

[is; 

[16 

[17; 



CLEO Collaboration, A. Freyberger et al, Phys. Rev. Lett. 76, 3065 (1996). 



D. Kaplan, e-print [hep-ex/9705002 



J. Liu and L. Wolfenstein, Nucl. Phys. B289, 1 (1987); G. C. Branco and M. N. Rebelo, 
Phys. Lett. B160, 117 (1985). 

J. L. Diaz-Cruz and G. Lopez Castro, Phys. Lett. B301, 405 (1993); J. L. Diaz-Cruz, 
J. J. Godina Nava and G. Lopez Castro, Phys. Rev. D51, 5263 (1995). 

C. W. Bernard, J. N. Labrenz and A. Soni, Phys. Rev. D49, 2536 (1994). 

J. D. Richman, Proc. of the 28th International Conference on High-Energy Physics, 
Eds. Z. Ajduk and A. K. Wroblewski, (World Scientific, 1997) p. 143. 

M. Wirbel, B. Stech and M. Bauer, Z. Phys. C29, 637 (1985). 

T. P. Cheng and M. Sher, Phys. Rev. D35, 3484 (1987); ihid D44, 461 (1991); A. 
Antaramian, L. J. Hall and A. Rasin, Phys. Rev. Lett. 69, 1871 (1992). 



8 



Channel 



Exp. BR 



upper bound 



DO 
DO 
DO 



e^e 



< 1.3 X 10 

< 4.2 X 10 

< 1.9 X 10 



-5 



a 



< 4.0 X 10-3 
a'''' < 2.3 X 10-3 
a'^^ < 4.9 X 10-3 



7rOe+e- 



DO 
DO- 
DO 
DO 

DO 

DO 
D+ 
D+ ■ 
D+ — > n'^ 
D+ ^ /v +//+// 



7rO/i+//- 
7rO/i±eT 



< 4.5 X 10 

< 1.8 X 10 

< 8.6 X 10 

< 1.1 X 10 

< 5.3 X 10 

< 1.0 X 10 

< 6.6 X 10 

< 1.8 X 10 

< 1.1 X 10 

< 5.9 X IQ- 



(7^^ < 4.2 X 10-2 
a"'' < 8.6 X 10-2 
cr'^^ < 5.8 X 10-2 
a^^ < 0.16 
ct'^^ < 0.38 
a^^ < 0.16 
a^^ < 2.2 X 10-2 
< 1.2 X 10-2 
cr'*^ < 2.9 X 10-2 
a'''' < 0.15 



DO. 
DO- 
DO- 
DO 

DO- 
DO 
D+- 

D+ - 



^ pOg+e- 



< 1.0 X 10 

< 2.3 X 10 

< 4.9 X 10 

< 1.8 X 10 

< 8.3 X 10- 

< 1.2 X 10 

< 5.6 X 10 

< 1.1 X 10 



a^^ < 0.35 
a'*'^ < 0.57 
a/"-'- < 0.25 
a^^ < 0.48 
a^'^' < 1.14 
ai"^ < 0.40 
a'^'* < 0.39 
n'''' < 0.96 



Table 1. Bounds on Yukawa couplings from FCNC and LFV D meson decays. 
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